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Wavepacket reconstruction via local dynamics in a parabolic lattice 
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We study the dynamics of a wavepacket in a potential formed by the sum of a periodic lattice 
and of a parabolic potential. The dynamics of the wavepacket is essentially a superposition of "local 
Bloch oscillations", whose frequency is proportional to the local slope of the parabolic potential. 
We show that the amplitude and the phase of the Fourier transform of a signal characterizing 
this dynamics contains information about the amplitude and the phase of the wavepacket at a given 
lattice site. Hence, complete reconstruction of the the wavepacket in the real space can be performed 
from the study of the dynamics of the system. 

PACS numbers: 03.65.-w, 03.75.-b, 32.80.Pj 



The reconstruction of the wavefunction of a quantum 
system in the real space from experimentally accessible 
measurements is seldom possible. Even when the spa- 
tial probabihty distribution is accessible, phase measure- 
ments in general require the use of a completely differ- 
ent, and often incompatible, technique. This is essen- 
tially due to the fact that a measurement of the posi- 
tion of a quantum system is usually accompanied by a 
vifavepacket reduction that erases the phase information. 
In a one-atom Magneto-optical trap (prepared in a re- 
productive way), repeated measurements of the fluores- 
cence light allow reconstruction of the atom probability 
distribution with a precision of the order of the wave- 
length of the detected light, whereas a measurement of 
wavepacket phases needs a specific atom-interferometry 
setup. 

In contrast, the quantum dynamics of a system is 
highly sensitive to both the amplitude and the phase of 
the wavepacket. The idea at the root of the present let- 
ter is to use information on the dynamics of the system 
to determine simultaneously both these quantities, with 
a spatial resolution at least comparable to that of the 
methods mentioned in the preceding paragraph. The dy- 
namics exploited here is related to the so-called "Bloch 
oscillation" (BO) observed in the motion of a quantum 
particle in a spatially periodic lattice to which a constant 
slope is added, forming a "tilted lattice" . In such a poten- 
tial, instead of climbing down the slope of the potential, 
the wavepacket performs bounded periodic oscillations in 
real and momentum spaces. The origin of this oscillatory 
motion can be retraced to quantum interference among 
pieces of the original wave packet that are transmitted 
and reflected by the successive potential barriers of the 
lattice, in much analogy with a chain of optical interfer- 
ometers with a linearly varying effective index The 
coherence of the wavepacket is thus a pre-rcquisitc to the 
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observation of the BO. The frequency of oscillation (the 
"Bloch frequency" ) corresponds to the energy shift from 
one lattice well to the next, and thus to the slope of the 
potential. This behavior, theoretically predicted in 1934 
has been observed experimentally in solid-state su- 
perlattices 0], with cold atoms §, with Bose- Einstein 
condensates l^] , and also with photons, in arrays of opti- 
cal waveguides |2|. A theoretical study of the BO in an 
frame close to that adopted here can be found in Ref. . 

Let us consider a quantum system, the "atom" , of mass 
M, placed in an one-dimensional potential formed by a 
spatially-periodic lattice of period d, to which a parabolic 
term ax'^/2 is added. With cold-atoms, this can be ob- 
tained by using a far-detuned standing wave which cre- 
ates a sinusoidal lattice ||, superposed to an in- 
dependent, focalized laser beam, producing a parabolic 
term (a laser-atom detuning much greater than the nat- 
ural width of the concerned electronic transition reduces 
spontaneous emission, insuring the system to be conser- 
vative). The Hamiltonian of the system is 

^--^^-^ocos(2..) + i^ (1) 

where we introduced dimensionless units in which dis- 
tances are measured in units of the lattice period d, en- 
ergy is measured in units of ad^ [or, equivalently, time 
is measured in units of h/{ad'^)], M* = Mad'^/h^ is a 
reduced mass, and h—1. Without loss of generality, we 
shall use a sinusoidal lattice. 

The energy spectrum of the Hamiltonian Eq. (^ is 
plotted in Fig. |l| and shows two distinct kinds of states. 
For energies below a threshold value e^, the correspond- 
ing eigenstates are localized in symmetric lattice wells —n 
and n (where the lattice site at n = is tied to the sym- 
metry axis of the parabola). For each possible value of 
the energy, one finds a symmetric and an anti-symmetric 
quasi-degenerated states, which are noted rcsp. (t>n{^) 
and (j)^{x). We suppose that the depth of the lattice is 
such that there are only one symmetric and one anti- 
symmetric state for each couple of lattice wells n and —n 
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FIG. 1: Energy spectrum of the Hamiltonian Eq. (|l|). 
Column (1) represents the energies of locahzed states, with 
a parabolic distribution. Column (2) displays the energies 
of the delocalized states. The inset displays the wavefunc- 
tion's modulus for a typical delocalized state (a), and a (anti- 
symmetric) localized state (b), compared to the parabolic lat- 
tice (c). Vo = 90. 



{Vq = 90). The corresponding eigenenergies are noted 
+ and + For energies above El, one finds 
both localized and delocalized states. At high energies, 
the latter are similar to the eigenstates of the harmonic 
oscillator and their eigenenergies approximately display 
a linear dependence on n (although this is not easily seen 
in Fig. I). 

We suppose that the atom wavepacket has been pre- 
pared in a superposition of localized states only. For 
high enough lattice amplitude Vq, one can neglect the 
energy shifts (5f and (5^ and consider as a new 
basis the left and right quasi-eigenstates (t>±n{x) = 
[(j)^{x) ± (/)^(a;)] /\/2, which are localized in the lattice 
site ±n with energy e„ — (up to an additive con- 

stant). Moreover, the localized eigenstates (pn{x) are al- 
most invariant under a translation by an integer number 
of steps of the lattice: 0„(x) « (f>m[x — {n — m)d]. This 
property is exactly verified for a tilted lattice, and is a 
good approximation in the present problem as long as the 
parabolic potential do not perturb too much the trans- 
lational invariance of the lattice wells over the spatial 
extension corresponding to the atom wavepacket (that 
is, as long as Vb » 1 in our units). 

At time t = the wavepacket (a;) is centered at no 
and has a spatial extension An: "fo = X)n '^"'^n(^)- 
time t one has: 



^{x,t) = Cn exp 
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For instance, the packet average position is 



exp 



im{n— — )t (3) 



{Xm = (4>n \x\ 4>n+m) IS almost independent of n). 

Eq. (|^) shows that, given the evolution of a quantum 
system, its coherences c„c*_„ can be obtained from the 
complex amplitude of the oscillations at the Bohr fre- 
quency LUnm- 



= m{n — m/2) 



(4) 



The parabolic component of the potential produces non- 
degenerate Bohr frequencies. In a tilted lattice cor- 
responding to a constant force F, the coherences be- 
tween neighbor sites c„c*_]^ correspond to an oscillation 
with the degenerate (n-independent) Bohr (or Bloch) fre- 
quency lob = Fd/h (ujb = F in reduced units). The 
parabolic potential, on the contrary, has a local slope 
F — X and the related "local" Bloch frequency depends 
on n. The local Bloch frequency (LBF) is by definition 
the Bohr frequency associated to neighbor sites at n and 
n — 1 (to = 1): 



(5) 



The presence of such a frequency in the dynamics is a sig- 
nature of the coherence c„c*_]^ that can be used to track 
the wavepacket amplitude and phase, as shown below. 

Experimentally, it is more convenient to work in the 
momentum-representation, because the spatial ampli- 
tude of the BO is in general very small whereas, in the 
momentum space, the amplitude of the oscillation is of 
the order of h/d (the extension of a Brillouin zone of the 
lattice), which is easy to resolve with available experi- 
mental methods. For cold atoms, for example, stimulated 
Raman transitions jl^ can be used to that end. A great 
deal of information can be obtained by just monitoring 
the evolution of the probability for a given value p of the 
momentum, P{p,t). We thus work in momentum space 
and introduce the eigenstates (pnip) in p— representation: 

ipn{p) ~ — [ ex-p{ipx)(f>no [x - {n - no)] dx 
V 27r J 

= exp [i(n - no)p] (p) (6) 

where we used the approximate translation invariance of 
the real-space eigenfunctions 0„ {x) . The evolution of the 
momentum probability-distribution is then: 

P{p,t) = |*(p,i)|2 = |vP„„(p)P X 



1 + exp(iTOp) 



from which one can determine the dynamics of the system 



(7) 
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FIG. 2: Fourier transform of the probability of finding the 
system in the zero-momentum state, Po{t). The envelope of 
both frequency folds well reproduces the envelope of |»I'o(2;)] 
(see inset). The arrows indicate the components correspond- 
ing to the coherences between neighbor-sites Qo{n) — c„c*_i 
(n = 20,21) and next-to-neighbor-sites Qo (21) — C2ic*9. 
Vo = 90, no = 21, and An = 7. 



Fig. g displays the Fourier transform of the quantity 
Po{t) = P{0,t) and has been obtained as follows. The 
time-dependent probability distribution Po{t) is calcu- 
lated by numerical integration of the Schrodinger equa- 
tion with the hamiltonian of Eq. (|l|) and with the ini- 
tial wavepacket 'i>o{x) shown in the inset of Fig. || (the 
phase of c„ is arbitrarily set to e"^'^^^). The spectrum of 
Po{t) presents components at the Bohr frequencies con- 
tained in Eq. (^), from which the amplitude and phase 
of the coherences c„c*_„ can be deduced. The inter- 
pretation of the spectrum may be complicated by the 
fact that Bohr frequencies corresponding to energy dif- 
ferences higher than next-neighbor wells may show up. 
These frequencies, correspondiiig to m > 1 in Eq. (|4), 
are higher than the LBF [Eq. (|^)] typically by a factor 
m. In the general case, the spectrum is thus composed 
of "folds" described by an integer value of m. However, 
a judicious choice of position of the wavepacket in the 
parabolic potential allow to separate one fold from the 
others. Frequency folds m — 1 and m — 2 are resolved 
if no + An/2 +1/2 < 2 {no- An/2 - 1) or, roughly. 
An <^ no- In producing Fig. |^, we managed to satisfy 
this condition. The first frequency fold is clearly seen as 
the low-frequency part of the spectrum associated to the 
coherences c„c*_]^. In Fig. ^, the shape of the wavepacket 
is reproduced by the envelopes of both the first and the 
second frequency folds. This is the case if the wavepacket 
is smooth enough at the scale of the lattice step that we 
can take |c„| ~ |c„_i|, and then directly identify the am- 
plitude of the n-component of the spectrum with |c„|. 

The complete reconstruction of the spatial wavefunc- 
tion from the c„ coefficients requires, in principle, the 
knowledge of the eigenstates of system. The detailed spa- 



tial shape of the eigenstates may depend on experimental 
parameters that may not be known with a sufficient pre- 
cision. With our system, the shape of the localized states 
is almost independent of n, as it is essentially determined 
by the lattice, the parabolic term giving a small correc- 
tion. As they are sharply localized, they can be con- 
sidered as delta-function-like, probing the local value of 
the wavepacket. The detailed shape of the eigenstates is 
thus not required for the reconstruction, at the a preci- 
sion corresponding to a lattice step. 

We can now proceed to the complete reconstruction 
of a wavepacket. We consider the first fold (m = 1) of 
the spectrum, and the zero momentum probability Pq (t) 
only. The other folds have essentially the same struc- 
ture; other values of the probed momentum may be used 
to extract complementary information or improve the ac- 
curacy of the reconstruction. We thus reduce Eq. (0) to 
the simpler form: 



Mt) 



1 



-lexp 



n--}t 



(8) 



One then sees that the Fourier component Qoin) = 
CnC^-i of Qo at frequency n— 1/2 directly gives the coher- 
ence between neighbor sites. If the wavepacket envelope 
is smooth, the amplitudes can be directly obtained by 
taking |c„p ~ \Qo{n)\. If the phase also varies smoothly, 
the phase shift from the site n — 1 to the site n is also 
directly obtained from the phase of Qo{n), up to an un- 
avoidable, but unimportant, global phase factor. 

A more precise determination of the coefficients c„ is 
obtained by using the second frequency fold which, as one 
can easily convince oneself, gives the coherences c„c* _2 = 
(5g(n), and the identity: 



Qo(n)Qo(n + l) 
Qo("+l) 



(9) 



The phase shifts can be straightforwardly obtained 
from the Qo{n). We display in Fig. the coherences 
c„ deduced from the spectrum of Fig. ^ The agreement 
with the amplitudes of the original wavepacket is excel- 
lent. The figure also shows real and imaginary parts of 
the complex amplitudes in very good agreement with the 
original wavepacket. 

We numerically verified that the above method can 
be applied to wavepackets of arbitrary shape. The re- 
constructed state then reproduces the wavepacket at the 
scale of the lattice step, all finer features are lost. A large 
number of delocalized states can appear in the decompo- 
sition of the initial state, but their contribution to the 
Bloch-oscillation spectrum is to add a wide background 
whose amplitude is much smaller than the peaks associ- 
ated to the coherence between localized states. 

Consider now the case in which there are many atoms 
simultaneously interacting with the potential. If the den- 
sity of atoms is low, so that there is in average one atom 
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FIG. 3: Reconstruction of the initial wavepacket. Real and 
imaginary parts of Cn are shown in (a), evidencing the very 
good agreement between phases and amplitudes of the origi- 
nal wavepacket components ( x ) and the reconstructed values 
(□). The method correctly reproduces the phase mismatch 
between neighbors (/> = 7r/4. Curve (b) shows explicitly the 
amplitudes versus n and compares these values with the ones 
of the original wavepacket. 



per site (or less), the present method reconstructs the 
atom distribution. If the atomic density is higher, the 
method measures the ensemble average of the atomic 
correlation (c„c*_„j). This quantity is very important 
because it determines the ability of the system to display 
coherent quantum dynamics [Q. Finally, note that the 
present method can also be applied to a tridimensional 
lattice: adding a parabolic potential in a given direction 
of the lattice will probe the wavefunction along this di- 
rection. 

In conclusion, we have presented a method allowing 
reconstruction of the spatial wavefunction of a quantum 
system using the Fourier spectrum of the Bloch oscil- 



lations in a periodic lattice superposed to a parabolic 
potential. Our method is quite general: it does not de- 
pend, e.g. on internal properties of the quantum particle 
being studied, or on the shape of the lattice wells. A 
parabolic shape is not necessary provided that (i) the 
Bohr (or Bloch) frequencies of the system are not de- 
generated and (ii) the approximate translational symme- 
try of the localized states is preserved. The method re- 
quires a potential shape presenting both a non-constant 
and non- vanishing slope in the region of interest: the lo- 
cal slope induces local Bloch oscillation and the spatial 
dependence, in turn, leads to the position-sensitivity of 
the method. In principle, the method is applicable to 
a variety of systems, including cold atoms in light po- 
tentials, Bose-Einstein condensates and electrons in an 
adequately constructed superlattice. Applied to a Bose- 
Einstein condensate, whose Bloch oscillations have been 
recently observed the present method may be used 
to probe nonlinear collective interactions. The method 
is also able to probe the spatial entanglement of two or 
more atoms, which can be of interest in the context of 
quantum information and quantum computing. We shall 
treat these interesting applications in the future. 
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